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BACKGROUND -1 

• GFSSP stands for Generalized Fluid System Simulation 
Program 

• It is a general-purpose computer program to compute 
pressure, temperature and flow distribution in flow 
network 

• It was primarily developed to analyze 

- Internal Flow Analysis of Turbopump 

- Transient Flow Analysis of Propulsion System 

• GFSSP development started in 1994 with an objective 
to provide a generalized and easy to use flow analysis 
tool 
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BACKGROUND -2 

DEVELOPMENT HISTORY 

• Version 1 .4 (Steady State) was released in 1 996 

• Version 2.01 (Thermodynamic Transient) was released 
in 1998 

• Version 3.0 (User Subroutine) was released in 1999 

• Graphical User Interface, VTASC was developed in 
2000 

• Selected for NASA Software of the Year Award in 2001 

• Version 4.0 (Fluid Transient and post-processing 
capability) is released in 2003 
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COURSE OUTLINE 

1 . Introduction & Overview 

2. Graphical User Interface 

3. Mathematical Formulation 

4. User Subroutine 

5. Pressurization, Waterhammer & 
Conjugate Heat Transfer 

6. Tutorials (Afternoon) 
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NETWORK FLOW OR NAVIER STOKES ANALYSIS - 1 
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NETWORK FLOW OR NAVIER STOKES ANALYSIS - 2 


Navier Stokes Analysis 

• Suitable for detailed flow 
analysis within a component 

• Requires fine grid resolution to 
accurately model transport 
processes 

• Used after after preliminary 
design 


Network Flow Analysis 

• Suitable for flow analysis of a 
system consisting of several 
components 

• Uses empirical laws of transport 
process 

• Used during preliminary design 
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NETWORK DEFINITION - 1 

GFSSP FLOW CIRCUIT 
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NETWORK DEFINITIONS - 2 


Network: 

O Boundary node 
□ Internal node 
-D— Branch 


At boundary nodes, all dependent 
variables must be specified 


m- 


-□ □- 

Series 




-m 



»| 5 ] 


Parallel 





-m 


Combination 


At internal nodes, all dependent variables 


must be guessed for steady flow and specified for transient 
flow. 
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NETWORK DEFINITIONS - 3 


UNITS AND SIGN CONVENTIONS 


Units 

- Length 

- Area 

- Pressure 

- Temperature 

- Mass injection 

- Heat Source 


External (input/output) Internal (inside gfssp) 

- inches - feet 

- inches 2 - feet 2 

- psia - psf 

- °F - °R 

- Ibm/sec - Ibm/sec 

- Btu/s OR Btu/lbm- Btu/s OR Btu/lbm 


• Sign Convention 

- Mass input to node = positive 

- Mass output from node = negative 

- Heat input to node = positive 

- Heat output from node = negative 
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DATA STRUCTURE 
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'i ljfi 


MATHEMATICAL FORMULATION - 1 

MATHEMATICAL CLOSURE - 1 
Principal Variables: 


Unknown Variables 

Available Equations to Solve 

1 . Pressure 

1. Mass Conservation Equation 

2. Flowrate 

2. Momentum Conservation Equation 

3. Temperature 

3. Energy Conservation Equation (First or Second Law 

of Thermodynamics) 

4. Specie Concentrations 

4. Conservation Equations for Mass Fraction of Species 

5. Mass 

5. Thermodynamic Equation of State 
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MATHEMATICAL FORMULATION - 2 


MATHEMATICAL CLOSURE -2 
Auxiliary Variables: 

Thermodynamic Properties & Flow Resistance Factor 

Unknown Variables Available Equations to Solve 

Density 

Specific Heats Equilibrium Thermodynamic Relations 

Viscosity [GASP, WASP & GASPAK Property Programs] 

Thermal Conductivity 

Flow Resistance Factor Empirical Relations 
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MATHEMATICAL FORMULATION - 3 

BOUNDARY CONDITIONS 


• Governing equations can generate an infinite number 
of solutions 

• A unique solution is obtained with a given set of 
boundary conditions 

• User provides the boundary conditions 
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MATHEMATICAL FORMULATION - 3 

A TYPICAL FLOW CIRCUIT 


Notes: 

1 ) Number of Internal Nodes = 7 

2) Number of Branches = 12 

3) Total Number of Equations = 7 x 4 + 12 = 40 

4) Number of Equations Solved by Newton 
Raphson Method = 7 + 12 = 19 

5) Number of Equations Solved by Successive 
Substitution Method = 3x7 = 21 


Atmosphere 
14.7 psia 



Helium 
151 psia 
70 °F 

66 


V 137 ) 

67 


r -. x (138) 
\25 ") 

Jl 

68 
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Graphical User 
Interface 
(VTASC) 

I Input Data 
j File 

• Creates Flow Circuit 

• Runs GFSSP 

• Displays results graphically 


PROGRAM STRUCTURE 


Solver & Property 
Module 


User Subroutines 
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GRAPHICAL USER INTERFACE - 1 


MODEL BUILDING 
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GRAPHICAL USER INTERFACE - 2 

MODEL RUNNING 
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GRAPHICAL USER INTERFACE - 3 

MODEL RESULTS 
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RESISTANCE & FLUID OPTIONS - 1 




1. Pipe Flow 


2. Flow Through a Restriction 


Thick Orifice 


L or — 



Square Reduction 



Thin Sharp Orifice 



Where: 

D, = Pipe Diameter 

D 2 = Orifice Throat Diameter 



3. Non-Circular Duct 


4. Pipe Flow with Entrance 
& Exit Losses 


5. Thin, Sharp Orifice 



Square Expansion 



Rotating Annular Duct 




6. Thick Orifice 


7. Square Reduction 


8. Square Expansion 


9. Rotating Annular 
Duct 
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RESISTANCE & FLUID OPTIONS -2 



Where: 

L = Duct Length 
o=Duct Rotational Velocity 
D = Duct Diameter 




B 



10. Rotating Radial Duct 


Labyrinth Seal 



11. Labyrinth Seal 


h .v. 



12. Face Seal 
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RESISTANCE & FLUID OPTIONS - 3 

GASP & WASP 


Index 

Fluid 

Index 

Fluid 

1 

HELIUM 

7 

ARGON 

2 

METHANE 

8 

CARBON DIOXIDE 

3 

NEON 

9 

FLUORINE 

4 

NITROGEN 

10 

HYDROGEN 

5 

CARBON MONOXIDE 

11 

WATER 

6 

OXYGEN 

12 

RP-1 
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RESISTANCE & FLUID OPTIONS - 4 


GASPAK 


Index 

Fluid 

Index 

Fluid 

1 

HELIUM 

18 

HYDROGEN SULFIDE 

2 

METHANE 

19 

KRYPTON 

3 

NEON 

20 

PROPANE 

4 

NITROGEN 

21 

XENON 

5 

CO 

22 

R-1 1 

6 

OXYGEN 

23 

R12 

7 

ARGON 

24 

R22 

8 

CM 

O 

O 

25 

R32 

9 

PARAHYDROGEN 

26 

R123 

10 

HYDROGEN 

27 

R124 

11 

WATER 

28 

R125 

12 

RP-1 

29 

R134A 

13 

ISOBUTANE 

30 

R152A 

14 

BUTANE 

31 

NITROGEN TRIFLUORIDE 

15 

DEUTERIUM ; 

32 

AMMONIA 

16 

ETHANE 

33 

IDEAL GAS 

17 

ETHYLENE 

34 

AIR 



35 

HYDROGEN PEROXIDE 
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ADDITIONAL OPTIONS 



• Variable Geometry Option 

• Variable Rotation Option 

• Variable Heat Addition Option 

• Turbopump Option 

• Heat Exchanger 

• Tank Pressurization 

• Control Valve 
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APPLICATIONS - 1 


FASTRAC TURBOPUMP 
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APPLICATIONS - 2 

GFSSP Model of the Fastrac Turbopump 


RP and He 
Drain 


XV 


Turbine 

Discharge 


RP Impeller Discharge 




O 01 


He Inlet 


> 


[202 | 




RP Impeller Inlet 
RP Inducer Inlet 


ftp* 

FS 


p** 


®FS 0 



604 

A 


P Pi 


240 


241 


2321 


2322 



P P 


2352 ■ 2351 



303 605 | 

A 

RP 

RP and Helium 
Helium 

02 and Helium 
02 


Resistance Types Used 

A Axial Annulus 
B Bearing 

FS Face Seal 

L Labyrinth Seal 

O Orifice 

P Pipe w/ Entrance and Exit 

Losses 


02 and He 
Drain 


02 Impeller Discharge 



02 Impeller Inlet 
02 Inducer Inlet 




Boundary Node 


R Rotating Branch 

T Branches Included Only for 

Thrust Calculations 


Internal Node 


Branch w/ Assumed 
Flow Direction 


* 6 pipes (5141, 5142, 5143, 5144, 5145, 5146) 

** 6 pipes (603 1,6033, 6034, 6035, 6036) 

*** 6 pipes (4181, 4182, 4183, 4184, 4185, 4186) 
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APPLICATIONS - 3 

Turbopump Test to 20000 RPM 
with Gas Generator 
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APPLICATIONS - 4 

Fastrac Turbopump Model Results 

Pressure history comparison at RP-1 Impeller back face [Labyrinth seal inlet (211) and outlet (212)] 



Time (s) 
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APPLICATIONS - 5 


LOX Tank 



RP-1 Tank 


Engine 

Interface 
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APPLICATIONS - 6 
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GFSSP Model of PTA Helium Pressurization System 


Facility Interface 
P=765 psia 
T=0-120 F 


Pipe SV7 Tee, Flo Thru Pipe SV3 Pipe Tee, Flo Thru Reduction Pipe SV4 

L=128 in. C l = 0.6 Kl=200 L=17in. C L =0.6 L=288 in. Kl=200 D^l.3 in. L=221 in. C L =0.6 

D=1.3 in. A=0. 63617 in 2 K2=0.1 D=1.3 in. A=0.63617 in 2 D=1.3 in. K2=0.1 D,=0.53 in. D=0.53 in. A=0.2827 in 2 


SV2 Pipe I Tee, Flo Thru 

C l = 0.6 L=19 in. If Kl=200 

A=0. 63617 in 2 D=1.3 in. K2=0.1 


-0 


v , .. Tee, Flo Thru 

Reduction —9 on 

D, =1.3 in. (| ^I 200 
D 2 =0.78 in. K2 01 




Reduction 
' 5 D,=l .3 in. 

" D,=0.53 in. 


A Pipe 
® I. 143 in. 
D=0.53 in. 


Flex Tube 
(&) L=28 in. 


65 

Engine Interface 
P=615 - 815 psia 


For stainless steel tubing, it is assumed 
that 8=0.0008 in. For steel flex tubing, 
this roughness is multiplied by four. 


Pipe 

L =Hin. (|) K2 =0.7 

D=0.78 in. 


SV15 
C l = 0.6 
A=0. 00001 in 2 ] 


Pipe A OF 15 
L=18 in. fc L =0.6 


(i035 y .- 

j SV5 
(§C L =0.6 
f A=0.63617 i 


(1037) Al (1041) 01 

Tee, Elbow Tee, Elbow 0 
Kl=500 Kl=500 J pi P e 

K2-0.7 I L=28 in - 


D=0.78 in. 


Pipe 
L=14 in. 
D=0.53 in. 


SV13 
C l = 0.6 
A=0. 00001 in 2 1 


SV14 
0 C l =0.6 

A=0. 63617 in 2 



f f. . . AoF13 OF 12 

J;_ 14 in -. {§C l =0.6 C l =0.6 

Jj- 0.53 in. A =0.00785 in 2 A=0.02895 ini 0=0.53 in. 


OF14 

C l =0.6 


Pipe 

I L-15 in. 



D=0.78 in. A=0.01767 in 2 A=0. 10179 in 2 D=0.78 in. 



Tee, Elbow 
Kl=500 FM2 

7 A C L =0. 77371 

K2_0 ' 7 A=0.55351 in 2 


Tee, Elbow 

(| C L -0. 83056 
K2 °' 7 I A=0.2255 in 2 

@ A0 0 


Expansion Pipe CV8 

Dj=0.53 in. L=14 in. C L =0.6 

D 2 =3.0 in. D=0.53 in. A=0.4185 in 2 


Expansion pipe 

D^O.78 in. L=21 in 


Propellant Tanks 
Control Valves 
Choked Orifices 
Various fittings 


CV9 

D=0.78 in. A=0.7854 in : 


Pipe Recovery 

r =0 6 * n- C L =0.0 ( 

Cl 0.6 j D=0.78 in. A=7.06858 in 2 


Recovery 
C L =0.0 
A=7. 06858 in 2 ] 


V ull =25 ft 3 


53 (1054) 54 

LOX Diffuser 
C l =0.6 

A=37.6991 in 2 


V prp =47 5 ft 3 

0L_ (1057). 58 

LOX Surface V A T 

C, -0.0 LOX Feedline 

A=4015. in 2 C L = dependent on 

stage of operation 
A=14.25 in 2 
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g«29> |0 |[jjj|- (103 

RP1 Diffuser RP1 Surface 

C l =0.6 C l =0.0 


A=37.6991 in 2 


A=3987. in : 


|0 LOX Pump 
Interface 


RP 1 F eedline ^ 

C L = dependent on 

stage of operation 
A= 14.25 in 2 0) 

RP1 Pump 
Interface 


Pressure (psia) 


Marshall Space Flight Center 
GFSSP Training Course 


APPLICATIONS - 7 

80 
70 
60 
50 
40 
30 
20 
10 
0 

-500 -400 -300 -200 -100 0 100 200 

Time (sec) 

GFSSP 5.0 Training Course 
Slide - 32 


Comparison of LOX Ullage Pressure with Test Data 


■ 




v v K 
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1 

1 
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L 1 \ ^ «, 



Jxf 

r» A A\ 1 
, 1 1 1 1 \ 1 
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1 
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Engine 
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Cut 
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APPLICATIONS - 8 

Verification of Fluid Transient Computation 


Liquid 
Oxygen at 
500 psia 
and 200 R 


Flowrate = 0.1 Ib/s 


D = 0.25 inch 


400 ft 


Valve closes in 100 
milliseconds 

450 psia 
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APPLICATIONS - 9 


Fluid Transient in Two phase flow 



Time(sec) 
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APPLICATIONS -10 


Predicted Fluid Transient Due to Condensation 
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SUMMARY - 1 

• GFSSP is a finite volume based Network Flow 
Analyzer 

• Flow circuit is resolved into a network consisting of 
nodes and branches 

• Mass, energy and specie conservation are solved at 
internal nodes. Momentum conservation is solved at 
branch 

• Generalized data structure allows generation of all 
types of flow network 

• Modular code structure allows to add new capabilities 
with ease 
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SUMMARY -2 

• Unique mathematical formulation allows effective coupling of 
thermodynamics and fluid mechanics 

• Numerical scheme is robust; adjustment of numerical control 
parameters is seldom necessary 

• Intuitive Graphical User Interface makes it easy to build, run and 
evaluate numerical models 

• GFSSP has been successfully applied in various applications 
that included 

- Incompressible & Compressible flows 

- Phase change (Boiling & Condensation) 

- Fluid Mixture 

- Thermodynamic transient (Pressurization & Blowdown) 

- Fluid Transient (Waterhammer) 

- Conjugate Heat Transfer 
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SUMMARY -3 



• GFSSP is available from NASA/MSFC’s 
Technology Transfer Office for US 
Government agencies and contractors 

• An Audio-Video Training Course is also 
available 

• More information about the code and its 
methodology is available at 
http://mi.msfc.nasa.gov/GFSSP/index.shtml 
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@ VTASC - AN INTERACTIVE 
PREPROCESSOR FOR GFSSP 


File Edit Advanced Help 

1 1 e* a B X 


nr 

m 

n 

\ 

3 

^0 H 

S \ / 

m u 2K 

yr / m 

n 7“"^ 0- — L S— 

f\ o " v 

Ul m Sc\ 

f n 

n'^ 39 

3 1 J 

(7, 4875,7.275) : Inch 
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BACKGROUND -1 

Visual Thermo-fluid dynamic Analyzer for Systems and 
Components (VTASC) is a program designed to 
efficiently build flow network models for use in the 
GFSSP program. 

- Visually Interactive 

• Eliminates pre-design of models 

• Immediate feedback on model 

- Self-Documenting 

• Hard copy of flow network 

• Bitmap image of flow network for inclusion into papers and 
presentations 
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BACKGROUND -2 

- Eliminates errors during model building process 

• Automatic node and branch numbering 

• Save and restore models at any point in the model building 
process 

• Robust 

- Pushbutton generation of GFSSP input file 

• Steady and Transient cases 

• Advanced features such as Turbopump, Tank Pressurization 
and Heat Exchangers 

- Run GFSSP directly from VTASC window 

• GFSSP Run Manager acts as VTASC/GFSSP interface 
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BACKGROUND -3 


- Post-processing capability allows quick study of results 

• Pushbutton access to GFSSP output file 

• Point and click access to output at each node and branch 

• Built-in plotting capability for transient cases 

• Capable of plotting through Winplot 

- Cross platform operation 

• Program written in C++ 

• Uses cross platform C++ GUI toolkit 
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CREATING A CHART IN WINPLOT -1 


i VTASC Version 2.03 - D:/GFSSP/2002 Code Mods/GFSSP4/2003 Class/tut5.vts 


File Edit Advanced Run Help 

II a & a ft GF55P 




r* 

m 

□ 

\ 


ur 

(6.7875,8.4875) : Inch 


□ H- 




►After completing a model run, select Winplot from the VTASC Run Menu 
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CREATING A CHART IN WINPLOT -2 



•From Winplot’s File Menu, select “Import ASCII” 
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CREATING A CHART IN WINPLOT -3 



Importing - winpltb.csv 


TIME, FI 2 ,F45 ,F56 ,F23 VI 2 ,V45 ,V5G ,V23 

SECONDS,LBM/S,LBM/S,LBM/S,LBM/S,FT /S,FT /S,FT /S,FT /S,LBF/ln2,LE 
0.1 00000E+00,0.1 00000E+01 ,0.1 00000E-01 ,0.1 00000E-01 ,0.1 00000E+01 ,0. 

0. 200000E +00,0. G7081 7E +00,0.1 281 4GE +03,0.1 28070E+03,0.GG5000E+00,i 
0. 300000E +00,0. 745686E -0G,0. 1 28068E+03,0.1 28070E+03,-.1 74455E+01 ,0. 
0.400000E+00,-.1 37788E -04,0.1 281 33E+03,0.1 28070E +03,0.1 03532E+01 , — 1 
0.500000E+00,-.1 388G8E -04,0.1 281 20E+03,0.1 28070E+03,-.G0871 2E-02,-.1 1 
0.G00000E+00,-.1 38G34E -04,0.1 28081 E +03,0.1 28070E+03,-.1 28870E+01 
0.700000E+00,-.431 785E -04,0.1 28081 E +03,0. 1 28070E +03,0. 349236E +00,-. ' 

0. 800000E +00,-. 433820E -04,0. 1 28081 E +03,0.1 28070E+03,0.300882E-0G,-.1 , 
0.900000E+00.-.434583E-04 .0.1 28081 E +03.0.1 28070E+03.-.375357E-08.-.1 -lJ 

||o Skipped rows 

jjl File index column 

fo File title row r- Permanent 

11 1 File 


Cancel 


Import 


•Use the Browse window to select the files you wish to import 
•The default GFSSP Winplot files are “winpltb.csv” & “winpltn.csv” 
•Selecting a file opens the Importing window. Click Import. 
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CREATING A CHART IN WINPLOT -4 



Ready | | |NUM 

•From the Parameter Selection window, select the data you wish to plot 
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DEMONSTRATION 
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VTASC DEMONSTRATION PROBLEMS -1 


Receiving Reservoir 



Pump valve 
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VTASC DEMONSTRATION PROBLEMS -2 


r 




Tank 

v = 10 ft 3 


Initial Conditions 

p. = 100 psia 
T. = 80°F 




d = 0.1 in 

~r~ 


Atmosphere 
p = 14.7 psia 
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VTASC DEMONSTRATION PROBLEMS -3 
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SUMMARY 



• VTASC is a flow network model builder for use with 
GFSSP 

• Flow networks can be designed and modified 
interactively using a “Point and Click” paradigm 

• Generates GFSSP Version 4.0 compatible input files 
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MATHEMATICAL FORMULATION 


Node 

izi 


Fluid 

Mixture 


Node 

( 

_i=± 

\ 


<5> 


Single 
Fluid 
k = 1 




\ ^ 

Node 

i 

t 

) * 

\ 

7 

< 


Fluid 

Mixture 


<5> 


Node 

_Lz_3 


Q 


rt|j = -rt]i 


Single 
Fluid 
k = 2 


Node 

i±4 


Alok Majumdar 

Propulsion System Department 

Marshall Space Flight Center 
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Content 



• Mathematical Closure 

• Governing Equations 

• Solution Procedure 
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MATHEMATICAL CLOSURE 


Receiving Reservoir 


Problem of a Steady State 
Flow Network 

• Given : Pressures and 
Temperatures at Boundary Nodes 

• Find : Pressures and 
Temperatures at Internal Nodes 
and Flowrates in Branches 


^ Pipe 
L = 1 500 ft 



D = 6 in. 

^ 8/D = 0.005 

0(WRT Gravity Vector) = 95.74 


Gravity Vector 


Gate 
Pump valve 


Primary Variables 

p 2 ,p 3 ,T 2 ,T 3 ,mn,m23,m34 

Secondary Variables 

Pi? P 3 "> Pi"* P3 


Pump 
A 0 = 30888 
B 0 = -8.067E-04 


Pipe 

L = 1500 ft 

Gate Valve D = 6 in. 

^ = 1000 s/D = 0.005 

= 0.1 0(WRT Gravity Vector) = 95.74 


12 


23 


- 34 y 



Supply Reservoir Pump Outlet 

P = 14.7 psia Valve Inlet 

T = 60° F 


Valve Outlet Receiving Reservoir 

Pipe Inlet P = 14.7 psia 

T = 60° F 


Legend 


Branch 


1 Boundary Node 

flc~| Internal Node 
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MATHEMATICAL CLOSURE 

Problem of an Unsteady Flow 

Network Primary Variables 


• Given : Pressures and Temperatures at 
Boundary Nodes and Initial Values at 

Internal Nodes Secondary Variables 


• Find : Pressures and Temperatures at 
Internal Nodes and Flowrates in Branches 
with Time. 


pV\ pip) 


/ \ 

Tank 

v = 10 ft 3 

Initial Conditions 

Pj = 100 psia 
T. = 80°F 


T 


Atmosohere 
p = 14.7 psia 
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MATHEMATICAL CLOSURE 
Principal Variables: 

Unknown Variable Available Equations to Solve 


1 . Pressure 

1 . Mass Conservation Equation 

2. Flowrate 

2. Momentum Conservation Equation 

3. Temperature 

3. Energy Conservation Equation 

4. Specie Concentrations 
(Mixture) 

4. Conservation Equations for Mass Fraction of Species 

5. Mass (Unsteady) 

5. Thermodynamic Equation of State 
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MATHEMATICAL CLOSURE 
Secondary Variables: 

Thermodynamic & Thermophysical Properties 

Unknown Variable Available Equations to Solve 




Density 

Specific Heats Equilibrium Thermodynamic Relations 

Viscosity [GASP, WASP & GASPAK Property Programs] 

Thermal Conductivity 

Flow Resistance 


Unknown Variable Available Equations to Solve 

1 . Friction Factor 1 . Empirical Relations 

2. Loss Coefficient 2. User Specified 
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GOVERNING EQUATIONS 



• Mass Conservation 

• Momentum Conservation 

• Energy Conservation 

• Fluid Species Conservation 

• Equation of State 

• Mixture Property 
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Coupling of Thermodynamics & Fluid Dynamics 
p - Pressure 

m - Flowrate 
h - Enthalpy 
c - Concentration 
p - Density 


Error 



Iteration Cycle 
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GOVERNING EQUATIONS 


MASS CONSERVATION EQUATION 


m 


t+At 


m 


At 


j = n 

= I m v 

j = 1 


Node 
j = 1 


Single 
Fluid 
k = 1 


Note : Pressure does not appear explicitly in 
Mass Conservation Equation although it is 
earmarked for calculating pressures 


Node 
j = 2 


Fluid 

Mixture 





\ ^ 

Node 

i 

( 

) > 

\ 


> 

< 


Single 
Fluid 
k = 2 


Node 
j = 4 


= 



Fluid 

Mixture 


Node 
j = 3 
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GOVERNING EQUATIONS 


MOMENTUM CONSERVATION EQUATION 

• Represents Newton’s Second Law of Motion 
Mass x Acceleration = Forces 


Unsteady 

Longitudinal Inertia 
Transverse Inertia 


Pressure 
Gravity 
Friction 
Centrifugal 
Shear Stress 
Moving Boundary 
Normal Stress 


External Force 
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MOMENTUM CONSERVATION EQUATION 


m, 


Mass x Acceleration Terms in GFSSP 
Unsteady 

(mu ) - (mu ) 

V ij / t+At \ ij / t 


U U 

o — — © — : — G> 


sA t 

Longitudinal Inertia 


MAX 


myfi 


(u.. - u ) - MAX - my ,0 ( u - u u ) 


o 
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MOMENTUM CONSERVATION EQUATION 


Force Terms in GFSSP 
Pressure 

ip, - p) A „ 


Gravity 

pgVCosO 

Friction 


-K mu 


Centrifugal 


m, 


A. 


pK rot CO A 
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GOVERNING EQUATIONS 

ENERGY CONSERVATION EQUATION 

• Energy Conservation Equation can 
be written in Enthalpy or Entropy 

• Based on Upwind Scheme 


Enthalpy Equation 

Rate of Increase of Internal Energy 


Single 
Fluid 
k = 1 


Enthalpy Inflow - Enthalpy Outflow + Heat Source 


m 


h - 


P_ 

pj 




f 


-m 


J 


h - 


r+Ar 


V 


p_ 

pj 


\ 


J 


At 


j=n 

= Z 

7=1 


Node 
j = 2 


Fluid 

Mixture 




Single 

Fluid 

k = 2 


Fluid 

Mixture 


Node 


Node 

> 

Node 

j = 1 



i 


j = 3 


tTlii = - 171;; 


< MAX 

1 

O 

€ 

1 

1 

h -MAX 

m ,0 

h > 

V 

IJ ' 

j 

ij ? 

i 


+ Q 
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GOVERNING EQUATIONS 

ENERGY CONSERVATION EQUATION 
Entropy Equation 


Node Fluid 
j = 2 Mixture 



Fluid 


Rate of Increase of Entropy = 

Entropy Inflow - Entropy Outflow + 
Entropy Generation + Entropy Source 


Mixture 



Single 

Fluid 

k = 2 


Node 
j = 4 


( /m )r + Ar Mr J MAX[- m t ,0 Js - MAX[m :j ,o]s 

Ar 7 — i " 7 " - 


J = n \ MAX\-m, g] l ' 

+ L i rvi r by gen 

7=1 1 kl J 


+ 


Q. 

i 

T 


i 


GFSSP 5.0 Training Course 
Slide - 65 


Marshall Space Flight Center 
GFSSP Training Course 


GOVERNING EQUATIONS 

FLUID SPECIE CONSERVATION EQUATION 


Node 
j = 2 

““ A 


Fluid 

Mixture 


Rate of Increase of Fluid Specie = 

Fluid Specie Inflow - Fluid Specie Outflow + 
Fluid Specie Source 


Single 
Fluid 
k = 1 


U,c M )„ sZ - 


A T 


j = n 

I 

7 = 1 



Single 

Fluid 

k = 2 


Fluid 

Mixture 


Node 


Node 

> 

Node 

j = 1 



i 


j = 3 


my = -nriji 


r 

< MAX 

1 

1 

o 

1 

c -MAX 

j,k 

1 

o 

1 

V 

J 

V 




> 


} + 


‘ i,k 
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GOVERNING EQUATIONS 

EQUATION OF STATE 

For unsteady flow, resident mass in a control volume is calculated 
from the equation of state for a real fluid 





R Tz 


Z is the compressibility factor determined from 
higher order equation of state 
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GOVERNING EQUATIONS 

EQUATION OF STATE 

• GFSSP uses two separate Thermodynamic Property Packages 
GASP/WASP and GASPAK 

• GASP/WASP uses modified Benedict, Webb & Rubin (BWR) 
Equation of State 

• GASPAK uses “standard reference” equation from 

• National Institute of Standards and Technology (NIST) 

• International Union of Pure & Applied Chemistry (IUPAC) 

• National Standard Reference Data Service of the USSR 
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GOVERNING EQUATIONS 


Mixture Property Relation 

Density 

• Calculated from Equation of State of Mixture with Compressibility Factor 


P = 


Pi 


k-n 


R = T.xR 


k = 1 


k k 


zikT , 

Compressibility Factor of Mixture is Mole average of Individual Components 

k-n 

Z. = I 


k - 1 


x z 
k k 


z k = 


Pi 


Pk^-kTk 


GFSSP 5.0 Training Course 
Slide - 69 


Marshall Space Flight Center 
GFSSP Training Course 


GOVERNING EQUATIONS 


Mixture Property Relation 


Thermophysical Properties 

• Viscosity, Specific Heat and Specific Heat Ratios are calculated 


by taking Molar Average 


k = n 

ft i ~ ^ %k ftk 
k = 1 


k-i i 

y i — 2 y fc 

k = 1 


k=n 



C „ X M, 
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GOVERNING EQUATIONS 

Mixture Property Relation 

Temperature 

• Mixture Temperature is calculated from Energy Conservation Equation 


£ Cp x* 7, MAX 

7=1 k=l 

r m , >°| 

+ (C mT) / At + Q 

\ p ,i i i / T ^ 1 

j=n k-n f 

EEC p lX ,MAX 

7=1 k=l 

m ,0 

- y 

1 + (C p .m) / At 


Limitation 


Cannot handle phase change of mixture 
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GOVERNING EQUATIONS 

Summary 

• Familiarity with GFSSP’s Governing Equations is not absolutely 
necessary to use the code 

• However, working knowledge about Governing Equations is 
helpful to implement various options in a complex flow network 

• A good understanding of Governing Equations is necessary to 
introduce new physics in the code 
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SOLUTION PROCEDURE 



• Successive Substitution 

• Newton-Raphson 

• Simultaneous Adjustment with Successive 
Substitution (SASS) 

• Convergence 
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SOLUTION PROCEDURE 



• Non linear Algebraic Equations are solved by 

- Successive Substitution 

- Newton-Raphson 

• GFSSP uses a Hybrid Method 

- SASS ( Simultaneous Adjustment with Successive 
Substitution) 

- This method is a combination of Successive Substitution and 
Newton-Raphson 
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SOLUTION PROCEDURE 

SUCCESSIVE SUBSTITUTION METHOD 

STEPS: 

1. Guess a solution for each variable in the system of equations 

2. Express each equation such that each variable is expressed in terms of 
other variables: e. g. X = f (Y,Z) and Y = f (X,Z) etc 

3. Solve for each variable 

4. Under-relax the variable, if necessary 

5. Repeat steps 1 through 4 until convergence 

ADVANTAGES: 

Simple to program; takes less computer memory 

DISADVANTAGES: 

It is difficult to make a decision in which order the equations must be 
solved to ensure convergence 
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SOLUTION PROCEDURE 

NEWTON-RAPHSON METHOD 

STEPS: 

1. Guess a solution for each variable in the system of equations 

2. Calculate the residuals of each equation 

3. Develop a set of correction equations for all variables 

4. Solve for the correction equations by Gaussian Elimination method 

5. Apply correction to each variable 

6. Iterate until the corrections become very small 

ADVANTAGES: 

No decision making process is involved to determine the order in which 
equations must be solved 

DISADVANTAGES: 

Requires more computer memory; difficult to program. 
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SOLUTION PROCEDURE 

SASS (Simultaneous Adjustment with Successive Substitution) Scheme 

• SASS is a combination of successive substitution and Newton-Raphson 
method 

• Mass conservation and flowrate equations are solved by Newton-Raphson 
method 

• Energy Conservation and concentration equations are solved by 
successive substitution method 

• Underlying principle for making such division: 

• Equations which have strong influences to other equations are solved 
by the Newton-Raphson method 

• Equations which have less influence to other are solved by the 
successive substitution method 

• This practice reduces code overhead while maintains superior 
convergence characteristics 




GFSSP 5.0 Training Course 
Slide - 77 


Marshall Space Flight Center 
GFSSP Training Course 



GFSSP Solution Scheme 


SASS : Simultaneous Adjustment 
with Successive Substitution 

Approach : Solve simultaneously 
when equations are strongly 
coupled and non-linear 

Advantage : Superior 
convergence characteristics with 
affordable computer memory 
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CONVERGENCE 



• Numerical solution can only be trusted when fully 
converged 

• GFSSP’s convergence criterion is based on 
difference in variable values between successive 
iterations. Normalized Residual Error is also 
monitored 

• GFSSP’s solution scheme has two options to control 
the iteration process 

- Simultaneous (SIMUL = TRUE) 

- Non-Simultaneous (SIMUL = FALSE) 
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CONVERGENCE 



Simultaneous Option 
• Single Iteration Loop 

- First solve mass, momentum and equation of state by the 
Newton-Raphson (NR) scheme 

- Next solve energy and specie conservation equation by 
Successive Substitution (SS) scheme 

- Solution is converged when the normalized maximum 
correction, A max is less than the convergence criterion 


Amax 


= MAX 



Z=1 


O/ 

<D 

z 


N e is the total number of equations solved by the Newton-Raphson scheme 
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CONVERGENCE 



Non-Simultaneous Option 

• Inner & Outer Iteration Loop 

- Mass, Momentum and Equation of state is solved in inner 
iteration loop by NR scheme 

- Energy and Specie conservation equations are solved in 
outer iteration loop by SS scheme 

- Convergence of NR scheme is determined/!^ 

- Convergence of SS scheme is determinedly 


Al, = MAX Ak, , A„ A, or A, 


A,, = MAX 


Zf 

i = 1 A/ 


etc. 
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Convergence Characteristics For Simultaneous Option 
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Comparison of Convergence Characteristics between Simultaneous and Non-Simultaneous Option 

in Converging-Diverging Nozzle 



0.00000001 
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SOLUTION PROCEDURE 

Summary 

• Simultaneous option is more efficient than Non-Simultaneous 
option 

• Non-Simultaneous option is recommended when Simultaneous 
option experiences numerical instability 

• Under-relaxation and good initial guess also help to overcome 
convergence problem 

• A lack of realism in problem specification can lead to 
convergence problem 

• Lack of realism includes: 

- Unrealistic geometry and/or boundary conditions 

- Attempt to calculate properties beyond operating range 
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USER SUBROUTINES 

Solver Module/function User Subroutine/function 



Alok Majumdar 

Propulsion System Department 

Marshall Space Flight Center 

alok.majumdar@msfc.nasa.gov 
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CONTENT 


• Motivation and Benefit 

• How they work 
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MOTIVATION AND BENEFIT 

• Motivation : To allow users to access 
GFSSP solver module to develop 
additional modeling capability 

• Benefit : GFSSP users can work 
independently without Developer’s active 
involvement 
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How do they work? 

- A series of subroutines are called from various locations of 
solver module 

- The subroutines do not have any code but includes the 
common block 

- The users can write FORTRAN code to develop any new 
physical model in any particular node or branch 

What users need to do? 

- Users need to compile a new file containing all user routines 
and link that with GFSSP to create a new executable 
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GFSSP PROCESS FLOW DIAGRAM 


Solver & Property 


Preprocessor Module 



User Subroutines 


New Physics 


• Time dependent 
process 

• non-linear boundary 
conditions 

• External source term 

• Customized output 
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DESCRIPTION OF USER 
SUBROUTINES 

Twelve User Subroutines were provided: 

- SORCEM: External Mass Source 


- SORCEF: 

External Force 

- SORCEQ: 

External Heat source 

- SORCEC: 

External Concentration source 

- KFUSER: 

New resistance option 

- PRPUSER: 

New fluid property 

- TSTEP: 

Variable time step during a transient run 
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DESCRIPTION OF USER 
SUBROUTINES 

- BNDUSER: Variable boundary condition during 


- USRINT: 

transient run (Alternative to history file) 

Provide initial values and steady state 
boundary conditions 

- PRNUSER: 

Additional print out or creation of 

additional file for post processing 

- FILNUM: 

Assign file numbers; users can define 
new file numbers 

- USRSET: 

User can supply all the necessary 

information by writing their own 


code 
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Solver Module/function 


User Subroutine/function 
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Solver Module/function 


User Subroutine/function 
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GFSSP INDEXING SYSTEM 


• Node and Branch Variables are stored in one- 
dimensional array 

• Node variables include: 

- Name 

- Pressure, Temperature, Concentration, Thermodynamic 
properties 

• Branch variables include: 

- Name 

- Flowrate, Velocity, Resistance coefficients, Reynolds 
number 

• Three subroutines are made available to Users for 
finding location and«di®e«nf©Fsa given node or 

^ Slide - 94 ^ 

branch 


Marshall Space Flight Center 
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NODE & BRANCH INDEX 

• User defined node names are stored in NODE-array. 

• NODE-array includes both internal and boundary 
nodes. 

• Total number of elements in NODE-array is NNODES 

• The internal nodes are stored in INODE-array. 

• There are NINT elements in INODE-array. 

• Branch names are stored in IBRANCH-array 

• There are NBR elements in IBRANCH-array 
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SUMMARY 


• User Subroutines can be used to add new capabilities that are 
not available to Users through Logical Options 

• New capabilities may include: 

- Introducing new type of resistance 

- Incorporating heat or mass transfer in any given node 

- Variable time step for a transient problem 

- Customized output 

• Checklist for User Subroutines 

- Identify subroutines that require modifications 

- Select GFSSP variables that require to be modified 

- Make use of GFSSP provided User variables in your coding 
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TANK PRESSURIZATION 


I 
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TANK PRESSURIZATION 


Predict the ullage 
conditions considering 
heat and mass transfer 
between the propellant 
and the tank wall 

Predict the propellant 
conditions leaving the 
tank 
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Pressurant 


V 



1 

Propellant to Engine 


Marshall Space Flight Center 
GFSSP Training Course 



TANK PRESSURIZATION 

ADDITIONAL PHYSICAL PROCESSES 


• Change in ullage and propellant volume. 

• Change in gravitational head in the tank. 

• Heat transfer from pressurant to propellant. 

• Heat transfer from pressurant to the tank 
wall. 

• Heat conduction between the pressurant 
exposed tank surface and the propellant 
exposed tank surface. 

• Mass transfer between the pressurant and 
propellant. 


Pressurant 


i 



Propellant to Engine 
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TANK PRESSURIZATION 

CALCULATION STEPS 


For each time step calculate 

• Ullage and Propellant Volumes 

• Tank Bottom Pressure 

• Heat Transfer between pressurant and 
propellant and pressurant and wall 

• Wall Temperature 

• Mass Transfer from propellant to ullage 


Pressurant 





Propellant to Engine 
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TANK PRESSURIZATION 

ADDITIONAL INPUT DATA FOR PRESSURIZATION 


Pressurant 


PRESS 

Logical Variable to Activate the Option 

NTANK 

Number of Tanks in the Circuit 

NODUL 

Ullage Node 

NODULB 

Pseudo Boundary Node at interface 

NODPRP 

Propellant Node 

IBRPRP 

Branch number connecting NODULB & NODPRP 

TNKAR 

Tank Surface Area in Ullage at Start, in 2 

TNKTH 

Tank Thickness, in 

TNKRHO 

Tank Density, lbm/ft 3 

TNKCP 

Tank Specific Heat, Btu/lbm - R 

TNKCON 

Tank Thermal Conductivity, Btu/ft-sec-R 

ARHC 

Propellant Surface Area, in 2 

FCTHC 

Multiplying Factor in Heat Transfer Coefficient 

TNKTM 

Initial Tank Temperature, 0 F 


4 



Propellant to Engine 


GFSSP 5.0 Training Course 
Slide - 101 


Marshall Space Flight Center 
GFSSP Training Course 


TANK PRESSURIZATION 

EXAMPLE 10 TANK SCHEMATIC AND GFSSP MODEL 


Helium 
95 psia 
120 deg. F 






Fluid: He 
P = 95 psia 
T = 120 °F 


1 



12 


C L = 0.6 
A = 0.785 in 2 


Aluminum Tank 
Diameter = 71.5 in. 

Wall thickness = 0.375 in. 


Ullage Node 


Pseudo Boundary Node 


Fluid: O 2 
P T =o = 74.76 psia 
T t =o = -264 °F 


Propellant Surface ; 34 ^ L = 


Propellant Node 


LOX Propellant to Engine 
P = 50 psia 


4 

A” 45 

5 


C T - 0.3043 



P = 50 psia 


A = 14.25 in 2 
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TANK PRESSURIZATION 

EXAMPLE 10 PRESSURIZATION INPUT 


NODE PRES (PSI) TEMP (DEGF) 

2 0 . 6700E+02 -0.2640E+03 

4 0.747 6E + 02 -0.2640E+03 

exlOhl . dat 
exl0h3 . dat 
ex!0h5 . dat 


MASS SOURC 
0 . 0000E+00 
0 . 0000E+00 


HEAT SOURC 
0 . 0000E+00 
0 . 0000E+00 


THRST AREA VOLUME CONCENTRATION 

0 . 0000E+00 0 . 4320E+05 1.0000 0.0000 

0 . 0000E + 00 0 . 8208E+06 0.0000 1.0000 


NUMBER OF TANKS IN THE CIRCUIT 
1 

NODUL NODULB NODPRP IBRPRP TNKAR TNKTH TNKRHOTNKC P TNKCON ARHC FCTHC TNKTM 

2 3 4 34 6431.91 0.375 170.00 0.20 0.0362 4015.00 1.00 -264.00 

NODE DATA FILE 
FNODE . DAT 
BRANCH DATA FILE 
FBRANCH . DAT 


Tank Input Units 

VOLUME, in 3 
TNKAR, in 2 
TNKTH, in 
TNKRHO, lbm/ft 3 
TNKCP, Btu/lbm-R 
TNKCON, Btu/ft-s-R 
ARHC, in 2 
TNKTM, deg. F 
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TANK PRESSURIZATION 

EXAMPLE 10 PRESSURIZATION OUTPUT 


SOLUTION 
INTERNAL NODES 


NODE 

P (PSI) 

TF (F) 

Z 

RHO 

EM (LBM) 

CONC 






(LBM/FT A 3) 


HE 

02 

2 

0 . 9138E+02 

-0 . 1347E+03 

0 . 1006E+01 

0 . 1047E + 00 

0 . 5144E+01 

0 . 9690E+00 

0.0310 

4 

0 . 98 69E+02 

-0 .2640E+03 

0 . 2310E-01 

0 . 6514E+02 

0 . 2 937E+05 

0 . 0000E+00 

1 . 0000 


BRANCHES 


BRANCH 

KFACTOR 

DELP 

FLOW RATE 

VELOCITY 

REYN . NO . 

MACH NO. 

ENTROPY GEN. 

LOST WORK 

(LBF- 

S A 2/ (LBM- FT ) 

A 2 ) (PSI) 

(LBM/ SEC) 

(FT/SEC) 



BTU/ (R-SEC) 

LBF-FT/SEC 

12 

0 .238E+05 

0 . 362E+01 

0 . 148E+00 

0 . 445E+03 

0 . 156E+06 

0 . 129E+00 

0 . 281E-02 

0 . 127E+04 

34 

0 . 000E+00 

0 . 000E+00 

0 . 163E+03 

0 . 899E-01 

0 . 412E+06 

0 . 114E-03 

0 . 000E+00 

0 . 000E+00 

45 

0 .263E+00 

0 . 487E+02 

0 . 163E+03 

0 .253E+02 

0 . 690E+07 

0 . 323E-01 

0 . 115E+00 

0 . 176E+05 


NUMBER OF PRESSURIZATION SYSTEMS = 1 

NODUL NODPRP QULPRP QULWAL QCOND TNKTM VOLPROP 

2 4 1.9642 8.5069 0.0022 196.4447 450.8641 


VOLULG 

49.1359 


SOLUTION SATISFIED CONVERGENCE CRITERION OF 0.100E-02 IN 5 ITERATIONS 
TAU = 10.0000 ISTEP = 100 


Tank Output Units 

QULPROP, Btu/s 
QULWAL, Btu/s 
QCOND, Btu/s 
TNKTM, deg. R 
VOLPROP, ft 3 
VOLULG. ft 3 
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TANK PRESSURIZATION 

EXAMPLE 10 ULLAGE AND TANK BOTTOM PRESSURE HISTORY 
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TANK PRESSURIZATION 

EXAMPLE 1 0 ULLAGE AND TANK WALL TEMPERATURE HISTORY 
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Mass Flow Rate (Ibm/s) 
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TANK PRESSURIZATION 

EXAMPLE 10 HELIUM FLOW RATE HISTORY 
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FLUID TRANSIENT 
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CONTENT 


• Classification of Unsteady Flow 

• Causes of Transient 

• Valve Closing 

- Comparison with Method of Characteristics 

• Valve Opening 

• Conclusions 
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CLASSIFICATION OF UNSTEADY FLOW 



• Quasi-steady flow is a type of unsteady flow when 
flow changes from one steady-state situation to 
another steady-state situation 

- Time dependant terms in conservation equation is not 
activated 

- Solution is time dependant because boundary condition is 
time dependant 

• Unsteady flow formulation has time dependant terms 
in all conservation equations 

- Time dependant term is a function of density, volume and 
variables at previous time step 
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CAUSES OF TRANSIENT 



• Changes in valve settings, accidental or planned 

• Starting or stopping of pumps 

• Changes in power demand of turbines 

• Action of reciprocating pumps 

• Changing elevation of reservoir 

• Waves in reservoir 

• Vibration of impellers or guide vanes in pumps or 
turbines 

• Unstable pump characteristics 

• Condensation 
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PROBLEM DESCRIPTION 



Rapid valve closing 



Objectives of Analysis: 

• Maximum Pressure 

• Frequency of Oscillation 
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Time (sec) 

Area (in 2 ) 

0.00 

0.0491 

0.02 

0.0164 

0.04 

0.0055 

0.06 

0.0018 

0.08 

0.0006 

0.10 

0.00 
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GFSSP Model 

L0 2 Propellant Tank 
500 psia, -260 ° F 


m 








• Run a steady state model with 450 psia ambient condition 

• Run unsteady model with steady state solution as initial value 

• Valve closing sequence is coded in BNDUSER 


9^0 

H 


10^1 



Valve 

Ambient 
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Description of Test cases 

Ptank = 500 psia 

Ttank = -260 0 F (Oxygen) = 70 ° F (Water) = -414 ° F (Hydrogen) 


Case 

No. 

Fluid 

Number of 
Branches 

Time Step 
(sec) 

Sound 

Speed 

(ft/sec) 

Flowrate 

(lb/sec) 

Pmax 

(psia) 

Period of 
Oscillation 
(sec) 

1 

lo 2 

10 

0.01 

2462 

0.0963 

626 

0.65 

2 

lo 2 

20 

0.005 

2462 

0.0963 

632 

0.65 

3 

lo 2 

5 

0.02 

2462 

0.0966 

620 

0.65 

4 

h 2 o 

10 

0.005 

4874 

0.071 

704 

0.33 

5 

lh 2 

10 

0.02 

3577 

0.0278 

545 

0.43 

6 

L0 2 & GHe 
(0.1%) 

10 

0.01 

1290” 

0.0963 

580 

1.24 

7 

L0 2 & GHe 
(0.5%) 

10 

0.01 

769” 

0.0963 

520 

2.08 

8* 

L0 2 (2 Phase) 
x exit = 0.017 

10 

0.01 


0.0963 

550 

1.17 

9’ 

L0 2 (2 Phase) 
x exit = 0.032 

10 

0.01 


0.0963 

538 

1.22 

10 

lo 2 


0.01 

2462 

0.0963 

611 

0.65 




Pressure oscillations are due to condensation 
Estimated from period of oscillation [a=4L/X] 


• Time step for each test case is so chosen that Courant number is less than 
unity 

• Courant Number = Length of Branch/(Sound Speed X Time Step) 


GFSSP 5.0 Training Course 
Slide - 114 


Marshall Space Flight Center 
GFSSP Training Course 


COMPARISON BETWEEN GFSSP & 
MOC SOLUTION FOR THREE FLUIDS 


Fluid 

Flowrate 

(lb/s) 

Velocity 

(ft/s) 

Friction 
Factor 
(Used in 
MOC 
solution) 

Sound 

Speed 

(ft/s) 

Max. Pressure 
rise above 
supply pressure 
(psi) 

Period of 
Oscillation 
(sec) 

MOC 

GFSSP 

MOC 

GFSSP 

Water 

0.071 

3.34 

0.0347 

4892 

214 

204 

0.33 

0.33 

Oxygen 

0.0963 

4.35 

0.0196 

2455 

136 

126 

0.65 

0.65 

Hydrogen 

0.0278 

19.01 

0.0157 

3725 

61 

45 

0.43 

0.43 
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COMPARISON BETWEEN GFSSP & MOC SOLUTION 



Z_l I I l_l I I I I I I I I I I I I I I I I I I I L 


12 3 4 

Time(sec) 
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PROBLEM DESCRIPTION 

Rapid Valve Opening 



Objectives of Analysis: 

• Maximum Pressure 

• Time to reach steady-state 


Time (sec) 

Area (in 2 ) 

0.00 

0.0 

0.01 

0.0088 

0.02 

0.1767 

0.03 

0.2651 

0.04 

0.3534 

0.05 

0.4418 
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GFSSP Model 


LH 2 Propellant Tank 
112 psia, -425 ° F 


HO 2 K— -0 -0 - 0 » 

12 23 34 45 

Valve 





Orifice 
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P (psia) 
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Pressure 



Time (sec) 
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T (R) 
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Temperature 
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Flowrate (Ibm/sec) 
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Flowrate 
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Quality 
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Quality 


1.00E+00 


9.00E-01 - 


8.00E-01 - 


7.00E-01 - 


6.00E-01 - 


5.00E-01 - 


4.00E-01 - 


3.00E-01 


2.00E-01 - 


1.00E-01 - 


O.OOE+OO 



O.OOE+OO 


5.00E-01 


1.00E+00 


1.50E+00 
Time (sec) 


2.00E+00 


2.50E+00 
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CONCLUSIONS 



• GFSSP has been used to compute fluid transient following rapid 
valve closure and opening 

• GFSSP predictions have been compared with MOC solution: 

- Maximum pressure and frequency compares well 

- Discrepancies exist in damping rate and shape of the pressure 
curve 

• Demonstrations have been made for 

- Sudden opening of cryogenic propellant in long pipeline 

• Time step must satisfy Courant condition 

• Predictions in all demonstration calculations show physical 
realism 
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Conjugate Heat Transfer 
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0 o 

1 I 

l l 


□ 


Boundary Node 
Internal Node 


Branches 


□ 

o 


Solid Node 
Ambient Node 
Conductor 


► 

<--> 


Solid to Solid 

Solid to Fluid 
Solid to Ambient 
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Mathematical Closure 



Unknown Variables 

1. Pressure 

2. Flowrate 

3. Fluid Temperature 

4. Solid Temperature 

5. Specie Concentrations 

6. Mass 


Available Equations to Solve 

1 . Mass Conservation Equation 

2. Momentum Conservation Equation 

3. Energy Conservation Equation of Fluid 

4. Energy Conservation Equation of Solid 

5. Conservation Equations for Mass Fraction of 
Species 

6. Thermodynamic Equation of State 
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Conservation 

Equation 
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SASS Solution Scheme 


Marshall Space Flight Center 
GFSSP Training Course 


Iteration Loop Governing Equations Variables 
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Verification of Conjugate Heat Transfer Results 
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Comparison with Analytical Solution 
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NBS Test Set-up of Cryogenic Transfer Line 
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GFSSP Model of Cryogenic Transfer Line 
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Comparison with Test Data 
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Summary 


• GFSSP has been extended to model conjugate heat transfer 

• Fluid Solid Network Elements include: 

- Fluid nodes and Flow Branches 

- Solid Nodes and Ambient Nodes 

- Conductors connecting Fluid-Solid, Solid-Solid and Solid-Ambient Nodes 

• Heat Conduction Equations are solved simultaneously with Fluid 
Conservation Equations for Mass, Momentum, Energy and 
Equation of State 

• The extended code was verified by comparing with analytical 
solution for simple conduction-convection problem 

• The code was applied to model 

• Pressurization of Cryogenic Tank 

• Freezing and Thawing of Metal 

• Chilldown of Cryogenic Transfer Line 

• Boil-off from Cryogenic Tank 
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